Unified parametrization of quark and lepton mixing matrices in tri-bimaximal pattern 
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Parametrization of the quark and lepton mixing matrices is the first attempt to understand the 
mixing of fermions. In this work, we parameterize the quark and lepton matrices with the help 
of quark-lepton complementarity (QLC) in a tri-bimaximal pattern of lepton mixing matrix. In 
this way, we combine the parametrization of the two matrices with each other. We apply this new 
parametrization to several physical quantities, and show its simplicity in the expression of, e.g., the 
Jarlskog parameter of CP violation. 



PACS numbers: 14.60.Pq, 12.15.Ff, 14.60.Lm 



I. INTRODUCTION 

Mixing of different generations of fermions is one of the 
most important unsolved problems in particle physics, 
and has already troubled physicists for half a century. 
Quark mixing of two generations was first introduced by 
Cabibbo [H decades ago as a unified explanation to the 
difference of the coupling constants in f3 decay and A de- 
cay in weak interaction. This idea was extended to the 
case of three generations by Kobayashi and Maskawa [3] ■ 
Parallelly, neutrino oscillation as the lepton mixing was 
originally motivated by Pontecorvo [3] fifty years ago to 
understand the solar neutrinoe deficit, and further gener- 
alized by Maki, Nakagawa and Sakata [11 . Visually, these 
mixings could be regarded as a rotation from fermion 
mass eigenstates to flavor eigenstates. Therefore, this 
rotation could be expressed as a mixing matrix, namely 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix V for 
quark sector and the Pontecorvo-Maki-Nakawaga-Sakata 
(PMNS) matrix U for lepton sector. 
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Both the CKM and PMNS matrices are unitary matrices, 
comprising three mixing angles 9ij (i, j = 1, 2, 3) and one 
Dirac CP-violating phase 6, and could thus be rewritten 
as 



C12C31 

-S12C23 - Ci2S23S3ie"' C12C23 - 
S12S23 - Ci2C23S3ie*'^ -C12S23 



S12C31 

Sl2S23S3ie 



iS 



iS 



S3ie 
S23C31 



diag(e*"i/^, e*"^/^, 1) multiplied to the matrix above. 
However, the two Majorana phases do not affect neutrino 
oscillations, and we ignore the diagonal matrix here. 

The similarity of the CKM and PMNS matrices arouses 
a natural question: what are their relations, e.g., the re- 
lations between their elements, or equivalently, mixing 
angles? Is there a unified comprehension of the two mix- 
ing matrices, just like the Cabibbo's way out to interpret 
different decays with the same coupling constant? Before 
the underlying grand unified theory of particle physics 
is discovered, phenomenological analyses play the main 
role in the correct direction. Amongst phenomenological 
analyses, to parameterize the CKM and PMNS matrices 
may be the first attempt. Parametrization provides us 
a numerical understanding of the hierarchy of different 
elements in the mixing matrix, and this is the purpose of 
our work. 

For the CKM matrix, various particle physics exper- 
iments have precisely fixed its value, with the allowed 
ranges of the magnitudes of its elements reading 
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We find that the CKM matrix is approximately the unit 
matrix with small deviations in the non-diagonal ele- 
ments, and it could be parameterized by the standard 
Wolfenstein parametrization Q, 
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where Sij = sin 



and Cij = cos 9ij . This expression 
was first introduced by Chau and Keung [s*] in 1984, 
and was advocated by the Particle Data Group in 1996. 
If neutrinos are of Majorana type, there should be one 
additional diagonal matrix with two Majorana phases 
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with A = 0.2272^ 



Nev- 



'T.'Z and A = 0.818iH?? a- 
ertheless, the situation of the PMNS matrix is exactly 
opposite. This is because (1) the long-existed inaccu- 
racy of neutrino experiments caused a delay in a practi- 
cal parametrization of the PMNS matrix as clear as the 
CKM matrix; (2) almost all the elements in the PMNS 
matrix are of order 0(1) except Ues, which makes the di- 
rect application of the Wolfenstein-like perturbative ex- 
pansion infeasible. 
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Whereas, circumstances have changed dramatically in 
recent years. Thanks to more and more precise mea- 
surements of neutrino oscillations, we have now a solid 
knowledge of the properties of neutrinos, aspects of their 
mass squared differences and mixing angles (although not 
as accurate as those of quarks) . This helps us to remove 
the first obstacle mentioned above, and in this paper, we 
focus on the second one, i.e., how to parameterize the 
PMNS matrix with large elements? Moreover, we try to 
find its relation with the CKM matrix with the help of 
the so-called "quark- lepton complementarity" (i.e., the 
sum of corresponding mixing angles in quark and lepton 
sectors is about 45°, always referred as QLC in litera- 
ture), and parameterize both of them. 

The organization of this paper is: in Sec. II, we take 
the tri-bimaximal pattern as the basis of the PMNS ma- 
trix, but do not parameterize the PMNS matrix around it 
immediately. However, we use the QLC to obtain a new 
basis for the CKM matrix, and discuss its parametriza- 
tion. In Sec. Ill, with the new parameterized CKM ma- 
trix, we utilize the QLC again, and the parametrization 
of the PMNS matrix around the tri-bimaximal pattern is 
consequently derived, with the same set of parameters in 
quark sector. Thus, we attain a unified parametrization 
of both the CKM and PMNS matrices. Several physical 
applications of the parametrization of the PMNS matrix 
are also studied in detail. After that, we talk about the 
relations between the two matrices and confirm the use- 
fulness of the QLC when we consider the parametrization 
of the PMNS matrix as an independent one. Conclusions 
and discussions are summarized in Sec. IV. 



parameters 


best fit 


2a 


3cr 


Amai 


7.6 


7.3 - 8.1 


7.1 - 8.3 


|Amii 


2.4 


2.1 - 2.7 


2.0 - 2.8 


sin^ 6*12 


0.32 


0.28 - 0.37 


0.26 - 0.40 


sin^ 6*23 


0.50 


0.38 - 0.63 


0.34 - 0.67 


sin^ 6*31 


0.007 


< 0.033 


< 0.050 



TABLE I: The three-flavor neutrino oscillation parameters 
from global data including solar, atmospheric, reactor (Kam- 
LAND and CHOOZ) and accelerator (K2K and MINOS) ex- 
periments. Amli is in units of 10~^ eV^ and |Am3i| of 
10-=^ eV^ 



(SS.Olj'j)" Atmospheric neutrino oscillation experi- 
mental data by different groups are: 1.5 x 10"'^ eV^ < 
lAm^gl < 3.4 X 



10-3 eV^ and sin^ 26*23 > 0.92 at 90% 



C.L. (Super-Kamiokande Collaboration ilflj); jAf 



^231 



cV^ and sin^ 2023 > 0.87 at 68% C.L. 



'23 



< 



2.74t°;44 X 10-3 

(MINOS Collaboration Tl^); 1.9 x IQ-^ eV^ < |Ar 
3.5 X 10-3 eV^ at 90% C.L. when sin^ 26*23 = 1, with the 
best fit value of 2.8 x 10"^ eV^ (K2K Collaboration [13]). 
On the other hand, the CHOOZ experiment [l3| indi- 
cated that the upper bound of the mixing angle ^31 is 
pretty small. (See Table U [l^ for a summary of the up- 
dated neutrino oscillation parameters and also Ref. [l5| 
for an excellent illustration of the oscillation parameters.) 
According to the results of global analysis, the modulus 
of the PMNS matrix is summarized as (at 3cr C.L.) [l^ 



II. PARAMETRIZATION OF THE CKM 
MATRIX 

Parametrization of the CKM matrix has already been 
well established, i.e., the Wolfenstein parametrization in 
Eq. (21), which explicitly show the deviations of the non- 
diagonal elements from the unit matrix at different pow- 
ers of the parameter A. While, this method was presented 
long before the release of enormous experimental data 
about neutrinos, and its parameters were introduced, in- 
dependently. Maybe a method of parametrization con- 
necting the quark sector to the lepton sector is more ef- 
fective at the epoch when fruitful neutrino experimental 
data are currently available. Therefore, in this section, 
we contemplate a new method for the parametrization of 
the CKM matrix, i.e., we parameterize it around a new 
basis, which is obtained from the tri-bimaximal pattern 
of the PMNS matrix by means of the QLC. Before con- 
tinuing this procedure, let us outline the main results of 
lepton mixing (i.e., neutrino oscillations). 

Nowadays, vast experimental data have confirmed 
the oscillations between different generations of neu- 
trinos. For solar neutrino oscillation, a global anal- 
ysis of the KamLAND experiments yielded Am^i = 
l.dtoi X 10-^ eV2 and tan^ 6*12 = OAOtom i; as for 
SNO experiments, Aml^ = S.OtoA x 10"^ cV^ and 6I12 = 



/ 0.77 - 0.86 0.50- 0.63 < 0.22 \ 
\U\^ 0.22- 0.56 0.44- 0.73 0.57- 0.80 , (3) 
\ 0.21 - 0.55 0.40- 0.71 0.59 - 0.82 / 

Eq. ([3]) manifests the nontriviality of the non-diagonal 
elements in the PMNS matrix, and illuminates a dis- 
tinct scheme for the parametrization of the PMNS ma- 
trix. We should not apply the method of the Wolfen- 
stein parametrization of the CKM matrix directly to the 
PMNS matrix, so firstly we need to find another basis 
for the PMNS matrix. Otherwise, the parameter used 
for the expansion would be of order 0(1), the expansion 
of the PMNS matrix then converges very slowly, and the 
hierarchy of different elements is also rather vague. 

The matrix chosen as the basis for our intention must 
be close to the experimental constraints in Eq. The 
most popular basis matrices are the "bimaximal" pattern 
(with sin 612 = sin 6*23 = \/2/2) and the "tri-bimaximal" 
pattern (with sin 6*12 = ^3/3 and sin 6*23 = \/2/2). The 
parametrization of the PMNS matrix based on the bi- 
maximal pattern was discussed in [l7|]. In this paper, we 
concentrate on the parametrization of the PMNS matrix 
around the tri-bimaximal pattern. 

Interestingly, the tri-bimaximal pattern was also intro- 
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duced by Wolfenstein [ig] thirty years ago with the form 



\/6/3 V3/3 
U= I -V6/6 V3/3 ^/2/2 
V6/6 -V3/3 V2/2 



(4) 



Recently, this special matrix has received dense dis- 
cussions in literature |19i] . Comparing Eq. ^ with 
Eq. ([4]), we find that the tri-bimaximal pattern agrees 
with the current experimental data fairly well, and the 
parametrization around it must be reasonable. 

Parametrization of the PMNS matrix around the tri- 
bimaximal pattern was treated in several approaches [20| . 
Here, we do not repeat these approaches, but postpone 
this parametrization in the next section. Now, we want 
to see what happens to the CKM matrix if we connect 
the tri-bimaximal pattern with the QLC. 

QLC is an appealing relation between the quark and 
lepton mixing angles. 



gCKM 



nPMNS 



^12 I ^12 
qCKM ^ flPMNS 
'31 ~ '^Sl 



= 45° 
« 0°. 



qCKM 

^23 



(5) 



Below we use the superscripts '^^^ and ^^^^ to denote 
the parameters of the CKM and PMNS matrices, respec- 
tively. This relation was first introduced by Smirnov 2l|, 
and has been heavily studied in different sorts of mod- 
els [12 . Eq. ([5]) is an extension for it. 

Now that we have these empirical equations, the CKM 
and PMNS matrices are not independent of each other, 
and therefore we are enlightened to parameterize them 
in the same framework. With these numerical equations 
and the tri-bimaximal pattern of the PMNS matrix, we 
easily get the trigonometric functions of the mixing an- 
gles in the CKM matrix, 

• /iCKM ~ 1 /iCKM ~^ ^ 

smt'in = ^— , cosOY-) = 1= — , 

sin0§r=O, cos0,T^ = l, 
sin0C™=O, cos03CK^ = l. 

The corresponding CKM matrix is thus fixed straightfor- 
wardly. 



/ 



V 




(6) 



Eq. ([6]), being the lowest order approximation of the 
CKM matrix, is not symmetric with 6*^2^^ = 9.7°, 
61CKM ^ 0° and e§^^ = 0°. Although this matrix is a 
little more complicated than the unit matrix, it is closer 
to reality. If we take the bimaximal pattern of the PMNS 
matrix with the QLC, we find that we get the unit matrix 
as the lowest order approximation of the CKM matrix. 
In the following, we just take Eq. ^ as the new basis 
for the CKM matrix, and parameterize the CKM matrix 
around it. 

Comparing Eq. ^ with Eq. ([1]), we make an expan- 
sion of V in powers of A (Attention, A here is not the 
Wolfenstein one, i.e., Vus = A), 



Vu. 



^/2- 1 

~7^~ 



+ A, 



where A — 0.0581^Q go|^g measures the strength of devia- 
tion of Vus from the corresponding element in the matrix 
we just obtained in Eq. ([6]). Since < A < 0.1, this 
expansion converges quickly. Similarly, we denote 



and A = 0.726lo;Qig. As the CP violating phase is in 
Kb, and Kb ~ 10 1 Kb I, we may set 



Kb = AX^{p-iTj). 



Thus, \/ p'^ + if = 1.6llo;gg, and the values of p and ij 
correspond to the Dirac CP-violating phase S'~^^^. 

In summary, four parameters: A, A, p and rj are intro- 
duced here, 



A = 0.058lt°:»l°, A = 0.726l°:°lt 



vVT^= 1.6ltO:|]5, 1.23 < 77/p < 4.70. 

where the ratio 77/p is fixed by rj/p ~ tan 7 = 

-tanarg[(KdK*b)/(KdK*b)] (7 = (63ti^)°) [6], and 
these four parameters characterize the CKM matrix com- 
pletely. 

Unitarity then determines the CKM matrix 



V = 



V2+1 V2-1 
V2-1 V2+1 




+ A 



-(30\/3- 21V6) 



-A2 



V2-1 



A^ 



-(3-2\/2) 1 

-1 -(3-2^/2) A 




V6 



2Ve 

^ip + ^v))A (. 



V6 



3_2^/2-^(p + ^?7) 



(30^/3 



V2+1 



A^ 
A 



(p - iri)A 




0(A3). 



(7) 
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We can extract the information of quark mixing from 
Eq. Q: 

1. The expansion is reasonable in powers of A and con- 
verges quickly. 

2. The term of A° is the new basis we have introduced 
in Eq. (O. 

3. The term of A^ clearly shows the deviation of the 
CKM matrix from this new basis, and it contains 
only one parameter A except A. 

4. The term of A^ is the modification of higher order, 
and the effect of CP violation exists in this order. 
Since CP violation is stored in Vub, the degree of 
CP violation in quark sector is of order A^ in our 
parametrization. 

In this section, the CKM matrix is parameterized in a 
new way, i.e., around a new matrix, which is considered 
as the basis of the CKM matrix. The crucial input is the 
QLC. Following this idea, we make use of the QLC once 
more, and the parametrization of the PMNS matrix is 
then obtained. 



III. PARAMETRIZATION OF THE PMNS 
MATRIX 

Although many properties of neutrinos have been 
known from neutrino oscillation experiments, informa- 
tion is not enough to ensure the accurate ranges of the 
modulus of the elements in the PMNS matrix. Known 
from Eq. ([3]), the parametrization of the PMNS matrix 
is troublesome if we do not take a basis matrix into ac- 
count. So we need a method to deal with it and make 
the parametrization feasible. Fortunately, the QLC con- 
nects the mixing angles of the two matrices which indi- 
cates that we may get the PMNS matrix from the CKM 
matrix. In this paper, we do not take the Wolfenstein 
parametrization of the CKM matrix as the perturbative 
expansion which leads to another parametrization of the 
PMNS matrix with the bimaximal pattern (see [23| for 
details), instead, we use Eq. ([7]) we have just obtained for 



the CKM matrix. With the QLC, we find the trigonomet- 
ric functions of the mixing angles in the PMNS matrix 
(to order of A-^), 



„• /iPMNS _ - " 



cos^P^^s 



sm 



cos 



gPMNS 
'23 



qPMNS 

'23 



V3 
3 
V6 
3 
V2 
2 

V2 
2 



(2%/2 - 2)A - (15%/6 - 21^/3)A^ 
(2 - V2)X - (15\/6 - 21\/3)A^ 



1-AX- -A^X^ 
2 



1 



AX - -A^X^ 
2 



sm 



iPMNS 
'31 

qPMNS 
'31 



Ax^ve+e, 

cost/Jf— " = (8) 
where A and A are the parameters we have defined for 
the CKM matrix, i.e., we describe the CKM and PMNS 
matrices with the same set of parameters. As we see in 
the parametrization of the CKM matrix, there are totally 
four parameters, but from Eq. ([5]) we only know two pre- 
cise numerical equations. It implies that we have to in- 
troduce another two new parameters ( and ^ to describe 
the PMNS matrix completely. From the fifth equation 



in Eq. 



sm 



iPMNS^-i5P™s 
^31 e 



where z = + i^. According to the experimental data 
of neutrino oscillations, the value of sin 6'™^^ has not 
been fixed yet, and what we have known is only its upper 
bound (see Table We find that \Ue3\ = sin6'ff < 
0.22 at 3ct C.L. from Eq. JS]). If 0.22 is taken for the value 
of \Ue3\, we find that \/ + ^ 90, which is unsuitable 
for the expansion of the PMNS matrix because the hier- 
archy of the first three orders would be indefinite. But 
if the value of \Ue3\ is smaller than 0.01, VC^TI^ < 4 
and sin6lfi'^NSg-»5i'"~s ^ j^^2^* preferable. The value 
of \Ue3\ depends on the experimental data, hence we 
may also set sin6lff NSe-»5™"^ ^ aX{C - i^) = AXz'* 
{^ '/(' ^ ^/^ = tanjP^^NS)^ ^ijgre z'* = C + ^C'- 
y/C'^ + ^'2 - 5 if we take 0.22 for the value of \Ue3\, 
which means that sin 0PMNS^-^s™^'^ ^ ^^^^f* .^^o^ld be 
better if the value of \Ue3\ is not small enough. Let us 
discuss the difference between the two cases. 



J 



Case 1: sin0™NSg. 



AX^z* 



With the trigonometric functions of the mixing angles in the PMNS matrix we have the parametrization of the 
PMNS matrix 



Ve V3 V2 



+X' 



2-V2 
2~-V2-^^A 



-(2-V2)-^A -{V2-I 
^(15\/6- 21^/3) 




-(15\/6- 2lV3) 



15V3-f + (2-^/2- ^)a+^^2 -(l5V3-f) + (V2-l-^)A-^A2 

V2- 1 



VeJ^^ 2V3^ 



-OiX^). 



z*A \ 



(9) 
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Case 2: smOl^''^ e-^'""""" = A\z'* . 

In this case, the last two equations in Eq. 



are replaced by (to order of 



sin 6' 



PMNS 
31 



AAv/(C'2+e'2), cosC 



PMNS 



Repeating the former process, we have 







\/6 

[/ = I _k k ^ I +A 

\/6 _ \/3 V2 
6 3 2 

/ -(15V6- 21%/3) 



a/2 



V2-f 



V%/6 73 



-(2V2- 2) 
^-1 + (7I-A 



A 



2 
V2 



+A2 



15^/3- +m^+pA2 



- (15V3 -^+mA- pA^ 15^3 - ^+nA + qA'^ 



+ 0(A3), (10) 



where 



= 2-V2- 



n = V2-l 



2z' 



V2' 



2V6 \/3' 



9 = 



2\/3 \/6' 



Now we see the meanings of the expansions of U in 
Eq. dH) and Eq. (fTU]): 

1. The PMNS matrix [/ is expanded in powers of A 
which makes the expansion converge quickly, the 
same as the CKM matrix V. So the parameter A 
we choose for the CKM matrix is quite well. 

2. The term of A° is the zeroth order approximation, 
we find that it is just the tri-bimaximal pattern of 
the PMNS matrix, i.e., the matrix we want to take 
as the basis of the PMNS matrix. So the method 
of getting the expansion of the PMNS matrix from 
the expansion of the CKM matrix with the QLC is 
feasible. 

3. The term of A^ shows the deviation of the PMNS 
matrix from the tri-bimaximal pattern. We find 
that the deviation is real and just relates with A 
and A in Eq. On the other hand, it is complex 
and relates with all four parameters in Eq. ^TU\i . 

4. The term of A^ is the modification of higher order, 
and it is much more complicated in Eq. (|10p than 
Eq. dH). 

5. CP violation is described by Ues in the PMNS ma- 
trix, and we find that the degree of CP violation in 
lepton sector is of order A^ in the parametrization 
of case 1 but of order A^ in that of case 2. Since 
the choice of case 1 or case 2 depends on the value 
of |/7e3|, it implies that the size of CP violation is 
determined by the magnitude of sin 6*™^^. 

Now let us apply the parametrization of the PMNS 
matrix to some physical observables. 



The amplitude for neutrinoless double beta decay 
{{)vj3(3) is proportional to the quantity @ 



Since the two Majorana CP-violating phases are related 
here, we must take them into account. 
In case 1 (to order of A^), 



l(w/3/3)] 



3 

o 



-(54-38V2)A^ 



We find that only mi and 7712 affect the value of | (t7t,^^)i |, 
7713 vanishes here to order of A^. Furthermore, there is 
only one parameter A which in our parametrization shows 
up in |(m/3/3)i|- 

In case 2 (to order of A^), 



1(^^/3)2! = 



{mp0)^+A^X^^{C'~iO^rn^ 



Apparently, \{mpf3)2\ is more complicated than |(rn^^)i|. 
Though the four parameters in our parametrization all 
appear in | (771^^)2 1, the effect of 7713 is still small, sup- 
pressed by A^. 

The rephasing-invariant Jarlskog parameter J [23| of 
CP violation is given by 

J = lin{Ue2U^3U*^U*2) = S12S23S31C12C23C31 sin (5. 
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An essential feature of J is that it is phase-convention 
independent which makes it very important when dis- 
cussing CP violation. The parameter is expressed in a 
simple form in our parametrization (to order of A^). In 
case we have 

J = ^AX^t 
6 



and in case 2, we have 




A AX^'. 



If we rescale (' AC and ^' — > A^, we will find that J' in 
case 2 becomes equivalent to in case 1 (to order of A^). 
If we take sin^ 0™^^ = 0.007 as a preliminary estimate, 
yJC'"^ + f 2 _ 2, then C' and are of order 0(1), case 2 



is preferable. However, if sin^ 6*™^^ = 0.0001 or less, C 
and ^ are of order 0(1), then case 1 is preferred. 

Leptonic unitarity triangles [1^ may be constructed 
with the orthogonality of different pairs of rows or 
columns of the PMNS matrix. The matrix elements C/e2, 
Ue5 and [7^3 play important roles in the experiments and 
they have relatively simple forms in our parametrization. 
There are two unitarity triangles which contain the three 
matrix elements, the 'V2, 1'a" triangle and the 'Vg, v^^^ tri- 
angle. Here we take the 'V2,i'3" unitarity triangle as an 
example (for previous discussions see Bjorken and King 
in Hi). 

The unitarity relation of the 'V2, 1^3" triangle is 



Ue2u:^ + u^2u;, + Ur2u;, = 0. 



(11) 



In case 1 (to order of A^), we have 

z 



Ue2U, 



e3 



u^2u;. 



V3 

1 

76 



AX' 



V2 



Ur2u;, 

and in case 2 (to order of A^), 



4_(1__)A+I^ 



Ue2U:, 



u,2u;. 



V3 

1 

V6 



AX 



2z' 



V2+I 
1 



V6 



AX^ 



V2 



45 
' 45 



63 



12 2V3 
63 z 



A 



V6 



A' X' 



V6 VT2 2\/3 



A 



V6 



A' X' 



1 



V2 2%/3 



.4 A 



45 



63 



112 



Ve VT2 V2 + 1 



A 



1 



Ur2u;, = - 1 - _ 



V6 



1 



V2 2V3 



63 



76 z' ^ \z'\ 
^ 3 ^/3 2^6^ 

^V6 z' , \zf 



A' \X' 



A+^--^ 



V3 



A' \X' 



We find that the Ue2Ul^ side of the 'V2,i'3" unitarity triangle is very small, suppressed by A^ in case 1 and A in 
case 2. This is because \Ue2\ is of 0(1), and |[/e3| is of 0{X^) in case 1 but 0(A) in case 2. Obviously, Eq. (fTT|) has 
been satisfied in both cases (to order of A^). If neutrinos are of Majorana type, there should be one additional factor 



multiplied to every expression, but Eq. pT|) still holds. 



oia2/2 

At last, we turn to neutrino oscillations. The probability of the oscillation from flavor Va to flavor P{va 
is given by 



Vc/*- 



''i 2E 



4=1 



Oa/3 



4 ^ Re{U*^^Up,U^,U;^) sin^ A,, + 2 ^ Iv^iW^^Up^U^.U;^) sin 2A, 



where a, j3 — e, ^, r, Aij = (m| — m'j)L/AE, L is the oscillation length and E is the beam energy of neutrinos. 

Let us take P{i'e v^) for example (for other oscillation probabilities see the Appendix). In case 1 (to order of 
A^), we obtain 



8 16 

9 ^ Y 



(\/6 - \/3)A - 8(\/2 - \)X^ ) sin^ A21, 
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and in case 2 (to order of A^), we have 

P(^, = 1 - (^^ - ^(V6 - V3)A - (^8(V2 - 1) + ^(C^ + e)A^) X') si 

-^(C" + r')^'A2(2 sin^ A31 + sin^ A32). 



sin^ A21 



r 



Certainly we find that P{va — > i^p) = 1 to order of 
A^, which means that the sum of the probabihty that 
a neutrino changes its flavor and the probabihty that it 
does not is unity. 

Nevertheless, in despite of the CKM matrix and 
the QLC, we may consider Eq. ([5]) as an independent 
parametrization of the PMNS matrix based on the tri- 
bimaximal pattern. Then Eq. ([9]) and Eq. pO]) are still 
obtained. According to the data in Table HI we find the 
ranges of the parameters in the PMNS matrix (at 2>a 
C.L.), 



-0.069 < A^'^^'^ < 0.079, 



-0.17 < A 



PMNS^PMNS 



< 0.16, 



and + C < 16 in case 1, + S,'^ < 1.3 in case 2. 



Without loss of generality, we set ^p^^^s ^ g^^^^ 
only acquire the upper bound of ^p'^ns^ 



A 



PMNS 



< 2.5. 



Clearly, the ranges of A^KM and ACKm 

are included in 

the ranges of A^^^^ and ylP^^^s gj^ce the experiments 
for the PMNS matrix are not as accurate as those for the 
CKM matrix. It implies that the corresponding param- 
eters of the two matrices may be correlated with each 
other to some extent when the parametrization of the 
PMNS matrix is independent. Let us see if there are any 
relations between the parameters of the CKM matrix and 
those of the PMNS matrix. 



sm 



qCKM I /)PMNS\ 
^12 + '^12 ) 



sin(0CK^+C3^NS) 



^/2 



(V6- V3)(A 



CKM 



- A 



-(9V2-12)AC™A™^^ 



V2 V2 
"2" ^ "2" 



PMNS 



) - (30 - 21 V2)(ACKM)2 _ (30^2 - 42)(A™^^)2 



^^CKM^CKM ^PMNS^PMNS^ 



V2 



CKM X CKM 



A 



PMNS ^PMNS^2 



It is clear that the QLC is reobtained to zeroth order, 
and the two equations above coincide with the first two 
equations in Eq. jS]) to order of A^ if A^^^ = A^^^^^ 
and ^CKM ^ ^PMNS^ If ^Yie relation between 9^^^ and 
^PMNS jg explicit, there should be an equation which links 
\l -V "(f' to a/C^ + ^2 , where the values of p and r\ cor- 
respond to b'^™^, and the values of C and ^ to (JP^^^s 
Thus, we state that the parameters in the PMNS matrix 
are not independent although the parametrization of the 
PMNS matrix is independent. The relations between the 
parameters of the CKM and PMNS matrices could be 
realized by the QLC. On another point of view, we may 
consider the QLC as the zeroth order approximation to 
the sum of the corresponding mixing angles of the CKM 
and PMNS matrices. 

The expansion of the PMNS matrix has been derived 
directly based on the tri-bimaximal pattern of the PMNS 
matrix in another work (see the first paper in [23j). In 
that case, the parametrization is independent, and the 
expansion is reasonable in powers of A (|A| < 0.1 at 3tT 
C.L.). The degree of CP violation is of order A^, the 



Jarlskog parameter is combined of two terms, and the 
effect of VTLz is suppressed by A^ in the amplitude of neu- 
trinoless double beta decay. It means that the expansion 
in that case is closer to case 2 in this paper where the 
parameters in the CKM and PMNS matrices are not in- 
dependent no matter the parametrization of the PMNS 
matrix is independent or not. The QLC plays an impor- 
tant role in the combination of the two matrices. Nev- 
ertheless, the magnitude of 0™^^ is not explicit, hence 
the parametrization of the PMNS matrix is arbitrary to 
some extent and the choice between Eq. ^ and Eq. pI7|) 
depends on the experimental data of |?7e3|. 



IV. SUMMARY 

In this work, we introduce a new basis for the CKM 
matrix, which is closer to reality than the unit matrix. 
The unit matrix is symmetric with no mixing angles while 
the new basis matrix is not symmetric with one non- 
vanishing mixing angle. Taking the tri-bimaximal pat- 
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tern as the zeroth approximation for the PMNS matrix, 
the new basis of the CKM matrix is obtained via the 
QLC. We parameterize the CKM matrix around the new 
basis instead of the unit matrix, and find the expansion 
is suitable as a perturbative approximation. 

The parametrization of the PMNS matrix is attained 
from the expansion of the CKM matrix with the QLC, 
then the two matrices are unified with the same set of 
parameters. Nevertheless, we get the expansion of the 
PMNS matrix, even if we consider parametrization in 
lepton sector as an independent one. In this case, how- 
ever, the parameters of the CKM matrix are still associ- 
ated with those of the PMNS matrix, since the QLC is 
reobtained when we combine the parameters of the two 
matrices. 

Anyway, the expansions of the CKM and PMNS ma- 
trices are reasonable and converge quickly since the pa- 
rameter A is very small. The efi'ect of CP violation is 
also seen from the parametrization of the two matrices. 



Furthermore, some observable quantities are expressed in 
relatively simple forms or only related with one or two 
parameters in our parametrization. Finally, the two ma- 
trices are correlated with each other and useful for ana- 
lyzing experimental data and phenomenology in particle 
physics. 
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APPENDIX 

In the appendix, we list all other probabilities for neutrino oscillations (to order of A^). 
In case 1, 



P{V^, Vr) 



^-^(V6-y3-A)A- 4(72-1) + ^ 



V6-V3-^Cl ^1 1 sin2 A21 



-^aA2(sin2 A31 - sin2 A32) + ^iAX^{sin2A^ 



\{V&-s/2 + A)\- \ A{s/2-l) 



1 — sin 2A31 + sin 2A32), 

2 ' sin^ A21 



l^fv6-^/3-fcV)A 



2V2 



CAA^(sin^ A31 - sin^ A32) 



-(Ve-TS- A)A 



72 




^^A2(sin2A2i — sin2A3i -I- sin2A32), 



V6- V3- 



V2 



C U A^ sin2 A21 



+—C,A>?{sa? A31 - sin^ A32) - — ^^A2(sin 2A21 - sin2A3i -h sin2A32), 



Piy^.-'V^) = 1- 



^-^(V6-y3-2^)A- (2(72-1) + ^ 



54- 38V2-t-^C^- 
54- 38V2+^CA+^A2 





\/3)A- f 2^2-2+^^2 



A^ sin2A2i 



+ Q - ^ ( V6 - ^)A + (^54 - 38^2 - \^ sin^ A31 



V3)A- (54 



38^/2 +^^2 



AMsin^A 



32 



\f2 

-h— ^AA^(sin2A2i - sin2A3i -|- sin2A32), 
o 



PK-i^e) = \^l-l{V6-V3 + A)X-(^4{V2-l)-^(v6-V3-y^CjAy'^^ 

/o /o 
^C^A^(sin2 _ gjjj2 ^^^^ _^ — ^AA^(sin2A2i - sin2A3i + sin2A32), 

Pii^r ^i^f.) = -(I'l^^- - - 2 + ^^'j A' j sin^ A21 



+ (I + ^(V6 - V3)A - (^54 - 38 V2 + ^A^^ X^^ sin^ A 



-^C^A^(sin2A2i - sin2A3i + sm2A32), 
PK-z..) = l-|^^-l(v^-V3 + 2^)A-(^2(v^-l)-^(^y6-V3-^c)A-^A2j^2j3i^2^^^ 

- (^^ - ^(%/6 - y3)A + (^54 - 38V2 - - ^A^^ A^ j sin^ A31 

- + ^(^6 - V3)A - [54 - 38V2 - + IaA xA sin^ A32. 



In case 2, 



- (^4{V2 ^l) + f(V6-V3)(^l- A + liC' + e)A'^ X'^ sin^ A,, 
+ (^^CAX - (^(^/6 - V3)CA - ^{C' + e)A'^ A^) sin^ A31 

- (^^CAX - (^(^/6 - V3)C'A + ^(C^ + e)A'^ A^j sin^ A32 
+ ^ (^'AA - (\/6 - \/3)^'^A2) (sin2A2i - sin2A3i + sin2A32), 



\{V2 - 1) - ^(V6 - V3) f 1 - ^C'l A + liC' + e')A'] X'] sin^ A21 



- (^^CAX - - ^/3)C'A + 1 (C^ + e)A'^ A^) sin^ A31 

- (^C'^A - (^(^/6 - x/3)C'A - ^(C^ + e)A'^ A^) sin^ A32 
~ {^'AX - (\/6 - \/3)^'^A') ( sin2A2i — sin2A3i +sin2A32), 
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- (^4{V2 - 1) + ^(V6 - x/3) (^1 - ^C') A + lie + e)A' 



^2 » A2 I sin^A2i 



^^"^C'AX -(^{V6- V3)CA + lie + e)AA sin^ A32 



3 ' \ 3 ' " 3~ 



3 



(sin2A2i — sin2A3i +sin2A32), 



9 9 \ V 
(^2iV2-l) + fiV6~V3) (l-^c) 



A-l(l + %'- le + MM AM sin' A21 



+ 1^54 - 38x/2 - ^(^^ - x/3)C'^ - ^ (l + ^/2C' - ^C" - J^") j j sin^ A31 



- 54 



_ 38^2 - ^(76 - V3)CA + l(^l-ylc + le + A'^ A^ j sin^ A32, 

i'(-.--r) = -(l-liV6-V3)x-Uv2-i) + l(i+^c+le + fe)AAx^]si^ 



3 3 



(^/6 - V3)A + (54 - 38x/2 - 1 (^1 + ^/2C' + ^C" + A') A^^ sin^ A, 

+ + ^ (^/6 - V3)A - (54 _ 38^ + ^ (1 _ |c' + ^C' + ^^'j 
+^ (e'^A - (\/6 - \/3)^'AA2) (sin 2A21 - sin 2A31 + sin 2A32), 



A2 aM sin^A 



32 



^+^(C''+C'VMAMsin2A2i 



3 

-(^4(^/2-l)-^(V6-V3) (l-^c) 
- (^C'AA- (^^(^-^)C'A+l(C'^+e'V'j A^j sin^A3i 
+ (^C'^A - (^(^/6 - ^/3)C'A - ^(C^ + e'^)^^ j A^ j sin^ A32 
+^ (e'^A - (^6 - V3)^'AA2) (sin 2A21 - sin 2A31 + sin 2A32), 
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sin^ A21 



1 _ 4 ^ V6 - 73) A + (54 - 38^2 - 1 (1 + x/2C' + \c'' + 1^ ]a']x') sin^ A31 



^ + ^(\/6-\/3)A- [ 54- 38^2+^ 



3 

'2 4 



(e'^A- (Ve- %/3)^'aa2) ( sin 2A21 — sin 2A31 + sin 2A32), 



(^54 - 38V2 + ^(^/6 - ^/3)C'^ - ^ (^1 + V2C - - Jc'^) 
[^54 - 38^/2 + ^(V6 - V3)C'A + ^ (^1 - ^C' + ^C'^ + 



A' sin^ A31 



'2 1 aM sin^A 



32- 
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